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Electronic energy transfer in photosynthesis occurs over a range of time scales and under a variety 
of intermolecular coupling conditions. Recent work has shown that electronic coupling between chro- 
mophores can lead to coherent oscillations in two-dimensional electronic spectroscopy measurements 
of pigment-protein complexes measured with femtosecond laser pulses. A persistent issue in the field 
is to reconcile the results of measurements performed using femtosecond laser pulses with physio- 
logical illumination conditions. Noisy-light spectroscopy can begin to address this question. In this 
work we present the theoretical analysis of incoherent two-dimensional electronic spectroscopy, 
2D ES. Simulations reveal diagonal peaks, cross peaks, and coherent oscillations similar to those ob- 
served in femtosecond two-dimensional electronic spectroscopy experiments. The results also expose 
fundamental differences between the femtosecond-pulse and noisy-light techniques; the differences 
lead to new challenges and new opportunities. 



I. INTRODUCTION 

In many optoelectronic devices, synthetic macro- 
molecules, and biological processes, electronic excitation 
energy is transported from one location to another. Elec- 
tronic energy transfer is thus an important process in 
physics, chemistry, biology, and engineering. A consider- 
able body of recent work has focused on electronic energy 
transfer in photosynthesis [1 j. The initial steps of photo- 
synthesis involve energy transfer over a range of distances 
and timescales: energy is transferred among the chro- 
mophores inside individual antenna complexes; from one 
antenna complex to another; and from the antenna com- 
plexes to the reaction center, where the excitation energy 
is used to separate charges and eventually split water. 
Two regimes for energy transfer in photosynthesis have 
been well understood since about the 1950s based on the 
balance between inter-chromophore electronic coupling 
(J) and the coupling of any individual chromophore to 
its bath (A) [2]. Forster theory can be used when the cou- 
pling to the bath dominates, J « A, and Redfield the- 
ory can be used when the coupling between chromophores 
dominates, J » A. But in the intermediate-coupling 
regime where J ~ A, energy-transfer processes are nei- 
ther well understood nor straightforward to treat math- 
ematically. Some theoretical progress has been made in 
the past two decades [3UTT]. 

Energy-transfer theories are important, and they need 
to be verified by experiments. Because the events are 
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fast — usually nanosecond timescales or faster — it is nat- 
ural for measurements to use femtosecond laser pulses. 
Most typical are pump-probe spectroscopy measurements 
[T2H26] , but technical advances [27-36 have made two- 
dimensional electronic spectroscopy (2D ES) measure- 
ments of photosynthetic proteins [37H46] frequent as well. 

One of the most interesting observations from 2D 
ES measurements of photosynthetic proteins are distinct 
cross peaks with oscillating amplitudes [37-40]. Cross 
peaks can be signatures of coupling (in several forms) 
or energy transfer [47J [48]. Cross peaks with oscillat- 
ing amplitudes are most often signatures of intramolecu- 
lar vibrational modes (vibrational coherence) or strong 
electronic coupling (electronic coherence) [4IJ [42j [491 
[52] . Similar coherent oscillations have been observed in 
photosynthetic proteins using related techniques as well 
[53H56] . The measured coherences have evoked ques- 
tions about energy transfer in the intermediate-coupling 
regime in biological systems and even larger questions 
about the nature of photosynthetic complexes [57H60] . 

Femtosecond spectroscopy is a powerful tool that will 
continue to be the primary source of new insight into the 
mechanisms governing energy transfer in photosynthesis. 
Yet a nagging question persists regarding the disconnect 
between the coherent excitation via femtosecond pulses 
used in experiments and the almost fully incoherent exci- 
tation via sunlight that photosynthetic organisms expe- 
rience in natural conditions: Can incoherent excitation 
produce the same coherences observed in femtosecond 
spectroscopy measurements? 

A few theoretical studies relate to this question [6T1 - 
165] . but no experiment has been performed or even sug- 
gested. Fortunately there is an established, joint theo- 
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retical and experimental framework that can address the 
question. This is the lesser-known technique called noisy- 
light spectroscopy which was developed concurrently with 
(one-dimensional) ultrafast laser spectroscopy. In 1984, 
Morita and Yajima [66], Asaka et al, [67], and Beach and 
Hartmann [68] independently demonstrated that noisy 
light could be used to achieve femtosecond-scale time 
resolution in degenerate four- wave mixing experiments. 
During the subsequent decade a variety of noisy-light 
analogues to more traditional ultrafast experiments were 
developed [6914101] . Although noisy light offers some 
advantages, it is not as generally useful as its short- 
pulse counterpart. One exception is coherent anti- Stokes 
Raman scattering spectroscopy (I^CARS), where the 
noisy-light version has proven to be a very useful tool 
that is still being used to investigate a variety of systems 
[87] I1Q2H116] . Noisy sources have found occasional uses 
and advantages in other spectral regimes as well [117- 
[121] . 

The main feature and utility of noisy-light spec- 
troscopy is that the time resolution is given by the co- 
herence time of the light, not the temporal envelope. In 
principle, the noisy beam may be continuous wave (cw), 
although in practice it is often pulsed on the ordered 
of nanoseconds (still essentially cw relative to femtosec- 
ond and picosecond material dynamics). In a typical 
experiment such as I^CARS, the noisy-light beam en- 
ters a Michelson interferometer to generate identical twin 
beams. One of the beams is delayed in time relative to 
the other by a controllable spatial delay in one arm of 
the interferometer. This can be generalized to multiple 
identical noisy beams with several controllable delays for 
more complicated techniques. 

Both short pulses and noisy light are spectrally broad, 
and in principle they could have identical optical spec- 
tral densities. But the differences between the excitation 
sources give rise to fundamentally different physical pro- 
cesses. For femtosecond pulses the frequency components 
must be phase locked (not to be confused with phase 
matched). That is, all frequency components must have 
a specific phase relationship so that the interference gen- 
erates a stable femtosecond laser pulse |122[ 1123] . Short 
pulses are ideally suited for direct time measurements; 
however, there are at least two disadvantages. First, 
femtosecond-pulse experimental setups are expensive and 
the experiments have strict stability requirements. In 
contrast, noisy-light experiments have more relaxed sta- 
bility requirements and require fewer expensive pieces of 
equipment. There are several reasons for this, but the 
most important is that noisy light is unaffected by dis- 
persion |112] , unlike short pulses where extensive efforts 
are often required [124] , Second, finely resolved spectral 
information must be determined through analysis of the 
time information. The phase-locking requirement pre- 
vents direct probing of the sample spectrally, although 
gross spectral probing is possible since experimentally ac- 
cessible femtosecond pulses do not have infinitely broad 
spectra. On the other hand, the spectrum of noisy light is 



completely phase unlocked. That is, the phase of any one 
frequency component is completely independent of any 
other frequency component. In a sense each frequency 
component is behaving as if it came from an indepen- 
dent cw source. In other words, the noisy light source is 
built from an incoherent superposition of monochromatic 
cw light; this random superposition produces a spatially 
coherent laser beam whose electric field is a stochastic 
function of time. 

It is worthwhile to consider how noisy light would in- 
teract with a photosynthetic pigment-protein complex. 
Because of its phase-incoherent nature, noisy light is an 
important step closer to sunlight than are femtosecond 
pulses. In this work we use the well-established tech- 
nique of factorized time correlation (FTC) diagram anal- 
ysis [110H112(I113[I125] — along with direct calculation of 
the spectroscopic signals — to describe a noisy-light ver- 
sion of two-dimensional electronic spectroscopy (I^ 4 ^ 2D 
ES, where I^ 4 ^ indicates four incoherent beams). The the- 
ory is presented and discussed for noisy light interacting 
with a Bloch four-level system. This system is complex 
enough to capture the essential physics of the interac- 
tion yet sufficiently simple to yield results that can be 
understood at a level giving significant physical insight. 

II. THEORY 

A. Input Fields 

The two-dimensional noisy-light experiment analyzed 
here involves three identical excitation beams, Ea, Eb, 
and Ec, arranged in the standard BOX geometry as indi- 
cated in Fig. []Ja). The fields produce a third-order non- 
linear response from the system under study. Beams Ea 
and Ec act in-phase while beam Eb acts out-of-phase; in 
the nomenclature of femtosecond 2D ES, beams Ea and 
Ec are 'nonconjugates' and beam Eb is the 'conjugate'. 
The resultant third-order signal is brought to quadrature 
at the detector with a fourth beam, the local oscillator, 
Elo- That is, the signal is heterodyne detected. As we 
detail below, the action of the LO is somewhat different 
than in femtosecond measurements. Working with the 
complex analytic signal, the temporal characteristics of 
the electric fields of the four input beams can be written 
as 

E A (t) = E p(t)e-^ (la) 
E B (t) = E p(t-a)e-^ t -^ (lb) 
E c (t) = E oP (t + r - a)e~ ) ( lc ) 
E L o(t) = E oP (t - ^e-*^-") , (Id) 

where Eo is a constant representing the field strength, 
p(t) is a complex stochastic function representing the 
random envelope of the noisy field and UJ is the central 
carrier wave frequency of the noisy light. The relative 
experimentally controlled temporal delays <r, r, £, and 
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FIG. 1. Parameters that would be used for I (4) 2D ES. (a) 
The BOX beam geometry. Noisy-light excitation fields Ea, 
Eb, and Ec generate a third-order nonlinear response in the 
sample, which emits a phase-matched beam in the direction 
of E s i g . The local oscillator field, Elo, is also in the direc- 
tion of i^sig- The four input noisy- light beams, Ea, Eb, Ec, 
and Elo, are identical except for the relative experimentally 
controlled time-delay variables a, r, £, and k. (b) The rel- 
ative pulse time orderings for femtosecond and 2D ES 
measurements, using the nonrephasing time ordering as an 
example. In femtosecond 2D ES, the experimental time-delay 
variables n , T2 , and T3 directly relate to the interaction times 
of the femtosecond pulses, £1, £2, and £3. In I^ 2D ES, the 
experimental time-delay variables have an ambiguous relation 
to the interaction times. The * indicates a common relative 
time point in all four noisy beams for ease of viewing. 
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FIG. 2. The energy-level scheme represents a pair of coupled 
two-level systems. We include all possible third-order signals 
produced by the transitions (/ii, [iy , /i2, and /x 2 ') among the 
four states: a common ground state, \g), two single-exciton 
states, \ei) and |e2), and a double-excitation state, |/). 

B. System 

In this work we calculate signals expected in a noisy- 
light spectroscopy experiment of the Bloch four-level sys- 
tem representing a pair of coupled two-level systems as 
indicated in Fig. [2] Coupling between the two-level sys- 
tems leads to an excitonic basis having energy shifts of 
the excited states and renormalization of the transition- 
dipole moments. In the exciton basis, there is a com- 
mon ground state, \g), two single-exciton states, |ei) and 
1 62), and a double-excitation state, |/), which can be 
reached via excited-state absorption from either of the 
single-exciton states. For completeness, we mention that 
this energy-level scheme could also represent a pair of 
uncoupled two- level systems. However, here we are in- 
terested in the coupled case, which means /ii 7^ \±y and 



k are shown in Fig. [T] Time-delay variable k is redun- 
dant, k = £ + <j — r, but it is introduced as a matter of 
convenience that will become clear below. We have also 
changed the sign of time-delay variable r to be opposite 
of what might be expected from femtosecond measure- 
ments to match previous noisy-light work and as another 
convenience that will become clear. We will also describe 
and show that although the light is on all the time, in the 
mathematical analysis we label specific interaction times 
and then integrate over them to account for all possi- 
ble interaction times. The relation between interaction 
times and the experimentally controlled time-delay vari- 
ables is simple for femtosecond 2D ES but complicated 
for l( 4 ) 2D ES. We depict this in Fig. [jjb), where experi- 
mental time-delay variables for the femtosecond measure- 
ment are given as n, T2, and T3 (in the literature there 
is an equivalent, alternative convention of r, T, and £, 
respectively) . 



C. Material Response 

The material response is handled via a Bloch model 
|126j such that 

i? (3) oc(-l)> (£)% (4) 

Xe -ifli(t 2 -ti) e -iCln(t3-t 2 ) e -itoiii(t-t 3 ) ^2) 

where £1, £2, and £3 are the times at which the fields 
interact with the material. The e -^(^- £ j) factors are 
the Liouville propagators during the intervals between 
field-matter interaction events, and Q x = uj x — ij x where 
uj x is the Bohr frequency of the coherence during the 
time interval between field interaction events, and 7^ is 
the phenomenological decay rate constant for that coher- 
ence. For population terms, Q xx = —iT xx1 where T xx is 
the excited-state lifetime of state x. A factor of (— l) n , 
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where n is the number of interactions on the ket side of 
the diagram, is also present to account for the sign dif- 
ference of the excited-state-absorption pathways. The 
factor of represents the product of four transition 
dipole moments specific to the Liouville pathway. 

For the four-level system depicted in Fig. [2j there are 
32 distinct Liouville paths that are triply resonant, 16 of 
which involve the higher energy state, /, and 16 which 
do not. Since noisy light is effectively cw, all orderings of 
fields Ea, Eb, and Ec must be considered. Because fields 
Ea and Ec are in-phase with one another, the 32 possible 
Liouville paths double to 64 when field ordering is taken 
into account. Each of the 64 Liouville pathways can be 
represented by a wave- mixing energy level (WMEL) dia- 
gram |102j or by a double-sided Feynman diagram [126 . 
Both sets of diagrams are shown in Fig. [3j Rephasing 
(Si), nonrephasing (5/j), and two-quantum (Sin) path- 
ways all potentially contribute to the signal at all delay 
times in the I^ 4 ^ 2D ES experiment, whereas in a fem- 
tosecond 2D ES experiment these sets of pathways are 
distinguishable by time ordering the fields. Because dif- 
ferentiation by time ordering is not particularly helpful, 
we have instead chosen to organize and label the groups 
of pathways according to the type of dynamics occurring 
during the second time period: G for ground-state bleach, 
P for population (stimulated emission), C for coherent 
oscillation, and T for two-quantum coherence. Excited- 
state absorption pathways carry an additional ' symbol. 

We then further differentiate the pathways using sub- 
script notation. The numeral subscripts I and II indi- 
cate Si and Sn pathways, respectively. Because all of 
the T pathways are by definition Sin pathways, they 
do not have a numeral subscript /// as that would be 
redundant. The numbered subscripts such as 1 and 2 in- 
dicate the first two or three transition dipoles following 
the convention set forth in Fig. [2] 



D. Detected Signal 

The third-order polarization is calculated in the same 
manner as for femtosecond nonlinear spectroscopy |126j 
by integrating over the fields and the material response 

Ef^cxP^cx f dt 3 f dt 2 f dhE A E* B E c R {3) . 

J— oo J —oo J— oo 

(3) 

In general all orderings of fields Ea, Eb, and Ec must 
be considered, although phase matching does restrict the 
possibilities somewhat. To maintain generality at this 
point the interaction time variables t\, t<i, and £3 are not 
associated with particular beams. To calculate a spe- 
cific term in the polarization, the time ordering of the 
fields must be specified. The third-order polarization is 
responsible for generating the signal field. The field is 
then brought to the intensity level through quadrature 



with the local oscillator field and averaged, 
I=(I(t)) = (\E LO + E^/) 

= (ho + ifi + E™El + Eif;E LO ) 
= (ho) + </<?>> + (Eif g Elo) + (c.c), (4) 

where I LO = E LO E* LO and if^ = E^E^, and we 
have suppressed the dependence of these expressions on 
the experimentally controllable time-delay variables. We 
consider that it is possible to remove term (Ig ig ) by ap- 
propriate balancing of beam intensities and to remove 
term {Ilo) through subtraction of a secondary measure- 
ment. Thus, suppressing the conjugate term, 

I^{E^ g El )^{P^El ) 

I oc ( [ dt 3 r dt 2 r dhEAE^EcEloR^) 

J — oo J — 00 J— 00 

cc [ dt 3 f 3 dt 2 t dt 1 (E A E* B E c E* LO )R^.(5) 

J — 00 J — 00 J — OO 

Other than the stochastic envelope terms and the as- 
sociated temporal averaging, the treatment thus far is 
similar to femtosecond methods; both are treatments of 
nonlinear optical signals in the semiclassical approxima- 
tion. 

Appropriate treatment of the averaging of the noisy- 
light process results in a four-point time correlator. Be- 
cause four-point time correlators are very difficult to 
handle mathematically, in noisy-light spectroscopy one 
assumes circular complex Gaussian statistics for the 
stochastic functions, p(t) |127j . This assumption — which 
has been used successfully in noisy-light spectroscopy for 
almost three decades — allows the four-point time correla- 
tor to be expressed in terms of two-point time correlators, 

(E A E B EcEl ) = (E A E B )(E c E LO ) 

+(E A E* LO )(E c E* B ). (6) 

This approximation makes the mathematics more 
tractable and still captures the essential features of noisy- 
light spectroscopy. Zhang and coworkers have considered 
noisy-light fields with a variety of statistical characteris- 
tics and have considered the four-point correlator directly 
for special cases [128 -130 . Normal treatment of higher- 
order correlators via the cumulant expansion leads to ex- 
ceedingly complicated analysis, which, to the best of our 
understanding, does not impact the experimental out- 
come. Example calculations beginning with Eqns. [5] and 
[6] are given in the Appendix. 

In short, the WMEL diagrams and the double-sided 
Feynman diagrams capture the light-matter interactions 
between the noisy fields and the Bloch four-level system 
that give rise to the third-order polarization that gener- 
ates the emitted signal field. This signal field is brought 
to quadrature at the detector with a local-oscillator field, 
which is also a noisy beam. Time averaging results in a 
four-point time correlator which reduces to two terms 
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FIG. 3. The sets of WMEL and double-sided Feynman diagrams for third-order noisy-light spectroscopy of the four-level 
system. The upper left, lower left, upper right, and lower right quadrants contain the stimulated-emission, ground-state bleach, 
excited-state absorption, and the two-quantum pathways, respectively. Shown are the pathways in which field Ea acts before 
field Ec, the a pathways. Not shown are the pathways in which field Ec acts before field Ea, the /3 pathways. The labelling 
scheme is described in the text. 



each having a product of a pair of two-point time cor- 
relators. As we will see below, each pair of two-point 
time correlators can be depicted and understood using 
another diagrammatic technique. 



III. RESULTS: FACTORIZED TIME 
CORRELATION DIAGRAM ANALYSIS 



The development of noisy-light spectroscopy has been 
hindered by the added mathematical complexity imposed 
by the quasi-cw and stochastic features of the light. The 
quasi-cw nature of the noisy light forces one to consider 
all possible time orderings of the field-matter interac- 
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tions in the perturbative treatment of the nonlinear sig- 
nal |131| . The stochastic nature of the light requires ex- 
plicit use of the bichromophoric model [104 to handle the 
nontrivial averaging at the intensity level. The total field 
is the phase-matched sum of the signal fields launched 
from each of the individual chromophores in the sample. 
In the homodyne detection scheme, the mod-square of 
this sum is dominated by the cross terms, meaning those 
fields launched from two distinct chromophores. The to- 
tal intensity is then very well represented by the sum 
of all pairwise, meaning two-chromophore, contributions. 
This bichromophoric model leads to the familiar N 2 de- 
pendence of the signal intensity in the coherent nonlinear 
spectroscopies, where TV is the number density of chro- 
mophores in the sample. For short-pulse spectroscopies, 
this point is not of practical relevance. The third-order 
polarization alone is sufficient for describing the nonlin- 
ear signal. Going to the intensity level is trivial — it is 
just the mod-square of the calculated polarization. For 
noisy-light spectroscopies, however, stochastic averaging 
of the noise at the intensity level requires that special 
attention be given to the signal intensity. The corre- 
lations among the various noisy-light field interactions 
on the two (otherwise independent) chromophores must 
be explicitly treated. These correlations are fundamen- 
tal to the understanding of noisy-light spectroscopies. 
Explicit use of the bichromophoric model requires as- 
signing distinct timelines (t and s) to each of the two 
generic chromophores. The two timelines allow for the 
(in general) different histories of evolution for each chro- 
mophore. Subsequent stochastic averaging links the t 
and s time variables in a nontrivial way. The heterodyne 
detection scheme considered for this work replaces the 
second chromophore field in the bichromophoric model 
with the local-oscillator field. In principle one would not 
need explicit use of the bichromophoric model but doing 
so is desirable because it mirrors the homodyne scheme 
and, as will be shown below, is beneficial for subsequent 
analysis. 



As it turns out, one can apply a diagrammatic rep- 
resentation of the mathematical expressions describing 
any given noisy-light signal |125j . These diagrams are 
called factorized time correlation (FTC) diagrams, and 
the complete set of these diagrams represents the break- 
down of the noisy-light signal into its elementary physical 
components. The set of FTC diagrams for the l( 4 ) 2D ES 
measurement considered here are presented in Fig. [4j In 
general, the set of FTC diagrams is isomorphic with the 
set of intensity-level terms for a given noisy-light spec- 
troscopy. Because the topological operations performed 
on the FTC diagrams represent mathematical operations 
performed on the analytic intensity- level terms, concep- 
tual tools have been developed to connect the topological 
operations with the mathematical operations. The con- 
ceptual tools also allow one to glean considerable physical 
insight from the FTC diagrams [TTOl [TT2l HT3l \T7E\ . 



A. Conceptual Tools 

Three of the conceptual tools — accumulation, synchro- 
nization, and color locking — will be particularly impor- 
tant for FTC diagram analysis of the I^ 4 ) 2D ES signal. 
Although extensive discussions have been previously pre- 
sented |112[ 1125] , these three concepts will be briefly re- 
viewed here because of their utility in the remainder of 
this work. 



1. Accumulation 

Across the ensemble of chromophores, the noisy fields 
may act on a given ensemble member at any time dur- 
ing the pulse duration. For those sequences of field ac- 
tions that produce the noisy-light signal considered in 
this work, the ability of the fields to act at any time on a 
given chromophore must be summed over the ensemble 
of chromophores contributing to the production of the 
nonlinear signal. From a topological point of view this is 
represented on an FTC diagram as the freedom of a tick 
mark to 'slide along a timeline'. Since the noisy beams 
are 'always' present, any time intervention (represented 
by a tick mark on the FTC diagram), or correlated pair 
of interventions, is free to slide along the timeline pro- 
vided the specific time ordering associated with a given 
FTC diagram is maintained. That is, during the sliding 
along the timeline, tick marks can not cross one another. 
To 'slide along the timeline' indicates the potential for 
the field action to take place at any time over which the 
tick mark is permitted to slide. However, an individual 
tick mark is locked to a partner tick mark by the segment 
representing a pair correlator. The two tick marks, thus 
linked, correspond to a correlated event pair that must 
slide along the timeline together. This ability to 'slide 
along the timeline' is called accumulation. Physically, 
accumulation represents the summation over the ensem- 
ble of chromophores in which the correlated noisy field 
interventions have the ability to act at any time relative 
to the other correlated pairs of actions that produce the 
signal. Mathematically, it represents integration over the 
time intervals between the light-matter interactions. 



2. Synchronization 

The second conceptual tool is concerned with the event 
coupling between the two timelines (t and s). This 
coupling involves inter-timeline pair correlators. This 
implies tight synchronization between the two chro- 
mophores of the events linked by a given correlator line. 
The 'precision' or 'strength' of this synchronization is in- 
versely proportional to the coherence time of the noisy 
light. For simplicity the correlation functions are taken to 
be (^-functions, thus synchronization is perfectly precise — 
the coherence time of the noisy light is taken to be zero. 
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F=F+ -^FF -#F-t- ^FF H^l— t- "^FF -F^=l — t- F=F} 



^ ^ ^ ^ ^ ^ ^ ^ 



T 'l2'2a T 'l2'2(3 T 21'la 1 21'1(3 1 21'2a 1 21'2(3 1 12'la 1 12'1(3 

-k-f^- -n^h +=cF^- -n^h -^d=^- -n^h -^d=^ 



T'i?'inf T'i 



FH=F F=F3- +^FF F=FF F^FF F=FF F^FF F^FF> 



Mri'a '12'1'P 1 21'l'a 1 21'1'P 1 21'2'a 1 21'2'P 1 12'2'a '12'2'f 

F^FF FH=F FFFF FH=4 FFFF +<-FF F=FF FF=% 



+H=F FFFF FrFF FFFF +H=F FFFF FH=F FFF^ 



FIG. 4. The set of 128 FTC diagrams. Although each FTC diagram represents a distinct mathematical term, many share the 
same or similar topological structure. The symbols <, o, and • at the ends of arrows represent time-delay variables a, r, £, 
and as, respectively, and the direction of the arrow indicates the sign of the time variable. 
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For this work synchronization occurs between the t time- 
line where the polarization on the chromophore is being 
developed and the s timeline representing the local oscil- 
lator field. 



3. Color Locking 

Noisy light has a broad spectrum similar to that of 
femtosecond pulses, with the important distinction that 
it is completely phase unlocked. One says the noisy light 
is color locked, because each color (or, equivalently, each 
frequency) is coherent only with itself — it is uncorrelated 
with any other color. Color locking is a consequence of 
the Wiener-Khintchine theorem |132j which is expressed 
mathematically most conveniently by examining a pair 
correlator in frequency space: 

{p{w)p*{o J '))=t{w)5{w-w') (7) 

where T(uS) is the optical spectral density of the noisy 
light. In the context of noisy-light spectroscopy, this im- 
plies that regardless of the spectral density of the sources, 
only identical frequencies may correlate to one another in 
a pairwise fashion. The consequence of color locking on 
FTC diagram analysis is extremely important: whichever 
frequency component happens to act from one field of a 
correlated pair (represented by one end of a line or ar- 
row), the other field of the correlated pair (represented 
by the other end of the line or arrow) must act with the 
same frequency component. This allows one to maintain 
the use of correlator terminology for a single color and 
refer to a x-x color-locked pair correlator, where x is a 
single frequency component of the noisy light. 



B. Heterodyne Detection 

To date, FTC diagram analysis has been applied to 
several noisy-light spectroscopies |113| [133-136 and also 
explored in a more general context |137( 1138] to elucidate 
the underlying mathematical structure of the diagrams. 
In all previous work the noisy-light signal was homodyne 
detected. This is the first FTC diagram analysis of a 
heterodyne-detected signal. The transition to the hetero- 
dyne case is accommodated in a straightforward manner. 
For the heterodyne case, the t timeline is drawn in the 
normal fashion as a horizontal line with three tick marks 
representing £ 1? and t%. (This is a condensation of the 
WMEL or double-sided Feynman diagram into a single 
line with tick marks.) The left hand side of the timeline 
represents — oo and the right hand side terminates with 
the quadrature event happening upon detection at time, 
t. The s timeline represents that of the local oscillator. In 
the heterodyne scheme the distinction between the two 
timelines is trivial but nonetheless formally important. 
There is a single tick mark placed at the right terminal 
side of the horizontal line representing the local oscillator 



action at the detector during the quadrature event which 
brings the signal to the intensity level. This means s = t. 

In femtosecond 2D ES, the local-oscillator delay is not 
related to the signal-generation process other than as a 
temporal reference that is removed via spectral interfer- 
ometry in post-processing of the data [139 . Importantly, 
for simulations of femtosecond 2D ES, the local oscilla- 
tor can be ignored altogether except mentioned that it 
allows one to recover the phase and amplitude of the 
emitted signal. In I^ 4 ^ 2D ES, the action of the local 
oscillator is slightly different. The local oscillator serves 
as a temporal anchor (setting s = t) because of synchro- 
nization (described above and in refs |112j and |125j ). It 
is interesting, but somewhat of an aside, that for homo- 
dyning in noisy-light experiments, the detection event 
is not part of the n-point field correlator. Hence the 
concept of a temporal anchor is not present in previous 
noisy-light works. Instead one builds up the conjugate 
signal (three field actions) separately on the s timeline. 
The n-point corrector (for the homodyne version of the 
measurement here it would a six-point correlator) and 
resultant pair correlators connect field events on and be- 
tween the t and s timelines. But these pairs are free 
to accumulate (slide along the timeline) and thus, al- 
though synchronized, do not form a temporal anchor on 
the laboratory time frame. Heterodyne detection thus 
has important implications both for the interpretation of 
the results and for experimental implementation as dis- 
cussed below. 



C. Analysis of the Topological Classes 

Inspection of the set of 128 FTC diagrams in Fig. [4] 
suggests a natural grouping based on topology. Three 
topological classes exist and are shown in Fig. [5] The 
line segments are labelled with the experimentally con- 
trollable delay variables: = ±a or ±r, (j) = ( or k. It 
is understood that for leftward pointing arrows 9 = — a 
and — r in all subsequent formulae. By way of exam- 
ple, a factor of the form e~ zQ0 would become e +2 ^ r for 
leftward pointing solid arrow and e~ zQ(7 for a rightward 
pointing open arrow. Generally, rephasing (nonrephas- 
ing) terms contain — (+0); all terms contain positive (j). 
It is interesting to note in Fig. [4] that a and £ only ap- 
pear together or r and k only appear together. There are 
48 unrestricted pathways, Fig. [5|a), 40 singly restricted 
pathways, Fig. |5|b), and 40 doubly restricted pathways, 
Fig.gc). 

One is able to perform FTC diagram analysis on each 
topological class to produce a single formula for generat- 
ing expressions for all diagrams in that particular topo- 
logical class. The set of 128 expressions can then be pro- 
duced with relative ease. In essence we can perform FTC 
analysis on three diagrams and obtain the analytic result 
for all 128 terms. For readers less familiar with FTC di- 
agram analysis, we reiterate that we have provided fully 
worked example calculations in the Appendix. FTC dia- 



t 2 e t 3 




FIG. 5. FTC diagrams capture the nontrivial relationship 
between the experimentally controlled time-delay variables (0 
and <fi) and the interaction-time variables (ti). Interaction 
times ti on the £ timeline are indicated, as well as the time 
delays and <j> as described in the text. The three topo- 
logical classes of FTC diagrams for this experiment are (a) 
unrestricted, (b) singly restricted, and (c) doubly restricted. 




FIG. 6. Analysis of the singly restricted FTC topological 
class. Because of the experimentally fixed time delays, the 
correlated pair of events at £2 and £3 are not able to slide 
along the timeline in an unrestricted manner — there is a lim- 
ited range over which they can slide. The directly correlated 
relationship between s and t\ and the directly correlated re- 
lationship between £2 and £3 result in an indirectly correlated 
relationship between £1 and £3 |112j . Indirect correlation is 
also present in this topological class between events £3 and s. 



in the unrestricted topological class are 



gram analysis is preferred because, in addition to reveal- 
ing the correct mathematical expression for each signal 
component, the analysis provides substantial physical in- 
sight and can help obviate mathematical errors. 



1. Unrestricted Topology 

We first consider the unrestricted topological class. For 
this class, Fig. [5|a), the line segments representing pair 
correlators are topologically disjoint. The intra-timeline 
segment is free to 'slide along the timeline' such that it 
can accumulate over the entire interval between £2 and 
£3. The inter-timeline segment synchronizes the last field 
event and the local oscillator. The intra-timeline seg- 
ment probes the response function on the interval be- 
tween the first and second field actions. This gives rise 
to a e~ %Vtld factor. The inter-timeline segment probes the 
interval between the third field action and the quadrature 
event with the local oscillator, giving rise to a e~ lVil11 ^ 
factor. The intra-timeline segment is free to accumulate 
fully over the entire response function between the second 
and third field actions. This contributes a factor 1/Qu 
to the signal. Assembled together — along with a factor 

of (|) (— l) n po/^ 4 ^o> where po is the 'dark' density op- 
erator factor, n is the number of ket-side interactions, 
M W is the particular sequence of four transition dipoles 
determined from the WMEL or double-sided Feynman 
diagram responsible for this term, and Iq is the cumula- 
tive intensity of the incident beams and local oscillator — 
the analytic expressions associated with FTC diagrams 



WM) = e[0]e[<A] 



(-i)Vo^ (4) / 2 



e — iQi9 e — iQui(f 



q 



(8) 

To obtain the term for a specific FTC diagram one needs 
only to read Qi, Qui, n, and p^ from the WMEL di- 
agram and and <fi from the FTC diagram. The results 
are collected in Table U 

Color locking plays an important role in this topologi- 
cal class. It forces the first and second field actions to be 
identical in frequency; therefore it is impossible to cre- 
ate coherent oscillations during the second time period 
for FTC diagrams in the unrestricted topological class. 
Thus all C terms in the unrestricted class are zero, as 
would be all T terms, however there are none. Only P 
and G terms persist. As we will show below, even though 
unrestricted terms cannot produce quantum beats from 
the C type terms, some of the unrestricted terms can 
produce oscillations during r due to polarization inter- 
ference. 



2. Singly Restricted Topology 

We consider now the singly restricted class. For this 
class, Fig. [5|b), the line segments are still topologi- 
cally disjoint. However, the inter-timeline segment syn- 
chronizes the first field action and the local oscillator. 
This has the effect of restricting the lower limit of ac- 
cumulation of the intra-timeline segment. The inter- 
timeline segment determines the range over which the 
intra-timeline can accumulate. Further, these FTC dia- 
grams represent nonzero terms only if <\> > \0\ because the 



10 



term 



GlUai 

Gmpn 

Gl\2ai 
Givifiii 
Gl21ai 
Gl21(3ii 
Gl22ai 
Gl22/3ii 
Plllai 

Pmpu 

Pl22ai 
Pl22f3ii 
P' 
P' 

Plllai 

Piiipu 



expression 



9[-a]9[C] 
9[-t]9[k] 
9[-a]9[C] 
9[-t]9[k] 
9[-a]9[C] 
9[-r]9[«] 
9[-a]9[C] 
9[-r]9[«] 
9[-a]9[C] 
9[-r]9[«] 
9[-a]9[C] 
9[-t]9[k] 
9[-a]9[C] 
9[-r]9[«] 
9[-a]9[C] 
9[-r]9[«] 
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n g e 1 <7 e -in e2 gt 


e+ 


^ge 1 t & — i^e 2 gK 


e+ 


Vgg 


e+ ' 


^gg 
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e - 




















n e2 e 2 
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^e 2 e2 


e + < 







term 
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GlH2ai 
Gui2(3ii 
GlI21ai 
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P' 

P' 

P'lIllOii 

P' 

Pllllai 
Pllllf3ii 

PlI22ai 
PlI22(3ii 



expression 



9[a]9[C] 
9[r]9[«] 
e[<r]9[C] " 
9[r]9[«] 
e[<r]9[C] 
9[r]9[«] 
9[a]9[<] 
9[r]9[«] 
9[a]9[C]^ 
9[r]9[«]^ 
9[a]9[C]^ 
9[r]9[«]^ 
9[a]9[C] 
9[r]9[«] 
9[a]9[C] 
9[r]9[«] 



~iCl ei g (<T + C) 
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in ei gT 9 ^- 


in & 2 g K 




^ ^gg 




e 




iQ-e-L g C 
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iCl e2 gr g J- 


i£l e -^ g K 



^gg 
in e2 g(cr+o 

^gg 
i ^e 2 g ( T + K ) 



l fe 2 " 



e iei 
neioT„-*"/ ei ' 



Je l e l 
e ig (<r + 







e 


i^l ei g (t + k) 






e - 


i ^ l e 2 g O + C) 




^ e 2 e 2 


e 




^ e 2 e 2 



TABLE I. Expressions for the 32 nonzero FTC diagrams with unrestricted topology. Each expression also gains a factor of 
(|) 3 ( — l) n pofi^ Iq . As detailed in the text, the 16 'C type terms with unrestricted topology are zero due to color locking. 
Interestingly, no terms in this class originate from two-quantum pathways. 



topology cannot accommodate <fi < \0\ (the diagram can- 
not be drawn). When <ft > the intra-timeline segment 
directly probes the response function between the second 
and third field actions via direct correlation. Indirect 
correlation |112[ 1113] is also present in these diagrams. 
The limited range of accumulation indirectly probes the 
response function between the first and second and be- 
tween the third and quadrature intervals. These FTC 
diagrams represent terms that offer the cleanest probe 
of coherent oscillations during time period r for 'C type 
terms. 

The FTC diagram analysis for this class shows that 
accumulation occurs over both the interval between ti 
and £2 and the interval between £3 and the quadrature 
event at s. As the intra-timeline segment slides along the 
timeline, the accumulation over each of those intervals is 



opposite from the other. This gives rise to a 



fac- 



tor. The interval between and £3 is directly probed by 
so one obtains a e~ lQlie factor. In addition to direct cor- 
relation, indirect correlation also appears as illustrated in 
Fig. [Hj This gives rise to a (er^^-d) _ e -ifi m (0-0)j 
factor. Taken together, this yields 



I S R = e[6}G[cj>-6} 



(-1)>oM (4) / 2 



-in ni {4>-e)\ 



Qui — Qj 



(9) 



This topological class is pathological because of the 
apparent pole when Vtm = fij. This is a removable pole, 



however, which is most easily seen by rearranging Eqn. 

mto 



i S R = e[e}G[0-e} 

x (l -e 



(-l)> /i (4) / 2 e-^e-^-e) 



in 



(10) 



which upon Taylor series expansion becomes, when 
Qui — >• and dropping higher-order terms, 



Isr = e[e]G[<t> - 8] 



(-1)>0M (4) ^ 2 e- m " 9 e- iQ '^ 



h 3 Qui — &i 

x(l-l + i(ft///-fi/)(0-0)) 



xe -iQ,(<j>-9) e -iQ u t 



(11) 



The mathematical analysis for such terms reproduces this 
result exactly, without the need for a Taylor expansion, 
because the factor of flm — fti appears in the argument 
of an exponential function before integration occurs. See 
the example in the Appendix for more details. The re- 
sults for this topological class are collected in Table [Til 



3. Doubly Restricted Topology 

Finally, we consider the doubly restricted topological 
class. For this class, Fig. [5|c), the line segments are not 
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TABLE II. Expressions for the 40 FTC diagrams with singly restricted topology. Each expression also gains a factor of 
(|) 3 ( — l) n pofi^ Iq . Interestingly, no terms in this class originate from rephasing pathways. 



topologically disjoint. The inter-timeline segment corre- 
lates the second field event with the local oscillator. Since 
the intra-timeline segment straddles the tick mark rep- 
resenting the second field event, the accumulation range 
is restricted from both the lower and upper limits. This 
topology is more complicated and must be piecewise de- 
fined about the point (/>=|#|.If|#|<(/> then the accumu- 
lation range of the intra-timeline segment is determined 
by its length, Fig. [7^a). Conversely, if |0| > </>, then the 
accumulation range is determined by the length of the 
inter-timeline segment, Fig. [7] (b). Another complica- 
tion is that the intra-timeline segment does not probe a 
single time interval exclusively but rather it probes both 
the interval between the first and second field actions 



and the interval between the second and third field ac- 
tions. These are indirect-correlation-based probes, so 
color locking does not apply [112 . Thus these FTC di- 
agrams can also probe coherent oscillations, albeit in a 
more complicated manner than those in the singly re- 
stricted class. 

FTC diagram analysis for both the \0\ > <p and the 
\0\ < (j) cases shows accumulation over all three intervals 

The indirect cor- 



which yields a factor of o — o ■ o — • 
relation probes of the response function differ for cases 
|0| > (j> and \0\ < (j) as illustrated in Fig. [7| For \0\ > (j), 
this gives rise to a factor of 
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TABLE III. Expressions for the 40 FTC diagrams with doubly restricted topology. Each expression also gains a factor of 
0[0]O[0] (|) 3 ( — l) n pQ Iq . Interestingly, no terms in this class originate from (one-quantum) nonrephasing pathways. 



and for \6\ < 0, this gives a factor of 
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Together, these factors yield 
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We used Heaviside functions to write the separate results 
as one large expression. See the example in the Appendix 
for further details. The results for this topological class 
are collected in Table Hill 
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FIG. 7. Analysis of the doubly restricted FTC topological 
class. Direct correlations, manifest as line segments, exist 
between t\ and £3 and between ti and s. Because the two 
segments are topologically linked, indirect correlation is com- 
plicated but nevertheless present. Most importantly, t\ and 
£2 are indirectly correlated, which relaxes the color- locking 
requirement. This enables coherent superpositions between 
the two single-exciton states to be developed on the 'C type 
terms. 



IV. DISCUSSION 



The expressions contained in Tables [l| [TTJ and III are 
the main result of this work. The total expected I^ 4 ) 
2D ES signal is the sum of these 128 expressions. In 
this section we present simulations for representative in- 
dividual terms and describe how the signals compare to 
both the femtosecond-pulse 2D spectrum [27H291 l32l l37l 
\M HH1 EH [T40l [HT] and the less-common frequency- 
domain 2D spectrum (a collection of transient-absorption 
measurements using a narrowband but tuneable pump), 
which has been performed in the infrared [142 -144 and 
ultraviolet |145j . but to the best of our knowledge not 
in the visible. We will relate our results to both the 
standard femtosecond 2D ES measurement — the sum of 
rephasing and nonrephasing components performed at a 
series of T2 times — and the two-quantum femtosecond 2D 
ES measurement |H1 [T46HT5T] . 

In this work, we use the Bloch four-level system to rep- 
resent a pair of coupled two-level systems characterized 
by the parameters (in the excitonic basis) given in Table 

IV We set = 1, and the normalized transition dipole 
moments were \i\ = 1, [12 = 1.2, jiy = 0.4, and ji2' = 0.5 
We computed the signals in the time domain using the 
given expressions; numerical Fourier transformation led 
to the presented I^ 4 ^ 2D ES. We describe several repre- 
sentative terms from each topological class that contain 
characteristics common to many terms. 

It is important to keep in mind that in principle all 
128 terms can contribute at all possible conditions of ex- 
perimental time delays. Thus we carefully consider the 
conditions under which each term is zero or nonzero. In 
multiple-beam spectroscopy measurements, positive and 
negative delay times can be easily confused. Our con- 
vention is depicted in Fig. [ljb) in which the field en- 
velopes can be either advanced or retarded by, for exam- 
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TABLE IV. Simulation parameters. Each frequency is defined 
as Q xy = 27r(h f xy — ij xy ) and its conjugate frequency defined 
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pie, retroreflectors on translation stages. The definition 
of positive a is such that path length is removed from 
field Ea so that a relative moment of its envelope inter- 
acts with the sample before that same relative moment 
from field Eb interacts with the sample. Time delay r is 
similarly defined on the interval between Eb and Ec-, ex- 
cept with the opposite convention; negative r means that 
a relative moment of the envelope of field Eb interacts 
with the sample before that same relative moment from 
field Ec- Finally, in the experimental implementation of 
the I^ 4 ) 2D ES measurement described here, time-delay 
variable £ must be scanned. This will probably be most 
easily accomplished by adding path length to the local 
oscillator beam; at each positive ( value where a rela- 
tive moment of the envelope of field Ec interacts with 
the sample before that same relative moment from field 
Elq, the intensity of the signal would be measured with, 
for example, a photodiode. Scanning ( is a necessary con- 
sequence that in noisy-light spectroscopy the envelope of 
the emitted signal has no deterministic time dependence. 



A. Unrestricted Examples 



In this topological class we study terms Gnu 
Giuai, Giiupu, Gi2iai, and G 



I21/3U- 



We begin with 



term Guuai, a term that originates from a nonrephas- 
ing Liouville pathway. We can write its expression as 



J(a,r,C) = e[a]e[C]^feP°e 
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to explicitly show the time-delay variables relevant to the 
measurement for this term. The Heaviside functions re- 
strict this term to be nonzero only for positive values of 
time-delay variables a and which means that this non- 
rephasing term is active under experimental time delays 
that correspond to a typical nonrephasing measurement. 
We present the time-domain signal in Fig. [8^a). 

Two-dimensional spectra are usually, but not always 
|152ffl57] , plotted as a function of frequency for the first 
and third time-delay variables. For the Gjmai expres- 
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In both the femtosecond and noisy-light measure- 
ments, the frequencies during the first and third time 
periods evolve phase of the same sign under exponential 
decay, which leads the peak for this term to appear in the 
(+, +) quadrant of Fig. |8jc). In addition, this term leads 
to a diagonal peak at the lower-energy exciton frequency 
in I^ 4 ) 2D ES just as the corresponding term would in 
femtosecond 2D ES. A significant difference between this 
component of the I^ 4 ^ 2D ES signal and the corresponding 
component of the femtosecond 2D ES signal is that this 
term does not depend on the waiting time, r, not even 
as a population decay. Instead, the population lifetime 
appears in the t^— pre-factor, which leads to the very 

1 99 

large amplitude of this peak. As will become clear, this 
is a common feature in I^ 4 ^ 2D ES. Because the signal is 
r-independent, this term is active even under conditions 
which seemingly violate causality: the probe sequence 
(fields Ec and Elo) can 'happen' before the pump se- 
quence (fields Ea and Eb) as long as the relative delay 
internal to each sequence is correct. This is a natural 
consequence of the fact that pair correlators — not indi- 
vidual field action times — govern dynamics in noisy-light 
spectroscopy measurements. 

The second main difference, as mentioned above, in- 
volves signal detection. In a direct experimental imple- 
mentation of I^ 4 ) 2D ES as described here, it will be nec- 
essary to scan the delay between fields Eq and Elo to 
extract in a point-by-point fashion the time dependence 
of the signal during the third time period. Implementing 
spectral interferometry in I^ 4 ^ 2D ES will be challenging 
because of the lack of causality and the asymmetry of 
the signal as a function of time-delay both of which 
are used in the spectral interferometry analysis in fem- 
tosecond 2D ES. 

The second unrestricted topology term we consider is 
Gnim, whose expression is 



FIG. 8. The result for term Gim^i. (a) The real part of the 
signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value I^ 2D ES spectrum. This term is 
independent of time-delay variable r, meaning its amplitude 
will not oscillate coherently nor will it decay with the excited- 
state lifetime. 



sion, Fourier transformation of a and ( yields 
PoM 4 /o 



where we inserted Q n 
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L 99 ~ ^99 ^ei9 ~ ^eig l leig- 

We present the absolute value (meaning amplitude) 2D 
spectrum for this term in Fig. [8]^c) on a linear color 
scale. There are several similarities — and two major 
differences — between the contribution from this term to 
femtosecond 2D ES and I^ 4 ) 2D ES measurements. 



I(a,T,()=e[- 



(15) 



and whose 2D spectrum is shown in Fig. [9j This term, 
which originates from a rephasing pathway, behaves in a 
manner similar to the femtosecond 2D ES measurement: 
it gains negative phase (u ge rather than uo eg ) as a func- 
tion of time-delay variable — a. This results in the peak 
appearing in the ( — , +) quadrant, which also happens in 
femtosecond 2D ES although most practitioners flip the 
sign of the excitation axis to positive frequency. This 
peak is also large in amplitude due to the t^- factor and 

1 99 

contributes to the lower-energy diagonal peak. Like the 
previous term, its amplitude is not r dependent. 

The third unrestricted topology term we consider, 
Gmipii, exposes another difference between the fem- 
tosecond and I^ 4 ) 2D ES experiments. This term has 
the same response function as G/nim> however, it is a 
f3 pathway, which accounts for the possibility that field 
Ec can interact before field Ea- In I^ 4 ^ 2D ES the con- 
tribution from this pathway must be considered because 
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FIG. 9. The result for term Gm a i. (top) The double-sided 
Feynman diagram, the WMEL diagram, and the FTC di- 
agram for this term, (bottom) Absolute value 2D ES 
spectrum. This term is also r independent. 



of the quasi-cw nature of the fields. The expression for 
this term is 



e[ T ]e[ K ]^jMA e -*n. 1B r e -tfi. 1Bf 
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where for the moment we suppress the time-delay vari- 
ables in the argument of the function. The time-delay 
variables natural to the expression are r and ft, but we 
must convert them into the variables we chose to rep- 
resent the total signal: a, r, and (. Using the identity 
r + n = <j + C, this change of variables results in 



7(a,r,C) = e[r]e[C + a-r] 



-iQ e 



19" e 



(17) 

This expression is almost identical to Eqn. [l3j however, 
the seemingly slight change in the arguments of the Heav- 
iside functions causes this term to behave in an unex- 
pected manner. Even though the term originates from a 
nonrephasing pathway, nonzero signal can accrue in the 
rephasing pulse time ordering — when a is negative — if r 
and ( are appropriately selected. 

The result is easiest to understand in the time domain, 
see the two plots in Fig. 10 Even though this term can 



be nonzero under rephasing conditions, it still accrues 
phase at positive frequency, therefore Fourier transfor- 
mation will produce a peak in the (+, +) quadrant. We 
display the expected contribution to the 2D spectrum for 
this term under both rephasing and nonrephasing mea- 
surements, when r = and when r = 0.2 ps, in Fig. 
The nonrephasing spectrum for r = is expected; 
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FIG. 10. The result for term Gimpa in the time domain 
(real part of the signal) when r = (left) and r = +0.2 ps 
(right). This term is zero when r < 0. The signal at r — 0.2 
ps is significantly weaker than the signal at r = 0. 
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the signal that appears under the rephasing conditions, 



FIG. 11. The result for term Gnwpa. (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (middle) Absolute value I^ 2D ES 
spectrum under nonrephasing conditions for two values of r. 
(bottom) Absolute value I^ 2D ES spectrum under rephasing 
conditions for two values of r. 



however, no has a very small anti-diagonal line width, 
and the maximum signal is almost an order of magnitude 
greater than the previously. The undesirable nonzero sig- 
nal under rephasing conditions can be eliminated in an 
experiment by proper Fourier filtering of the signal quad- 
rants. When r > 0, as in the 0.2 ps plots, the signal is 
also elongated about the diagonal but it is significantly 
weaker than the signal at r = under both rephasing 
and nonrephasing conditions. 

The fourth term we present is G//2im? which is dif- 
ferent from the two above because it involves transitions 
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FIG. 12. The result for term G II2 i a i- (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (bottom) Absolute value 2D ES 
spectrum. This term contributes with large amplitude to the 
cross peak and is independent of time-delay variable r. 



to both of the single-exciton states. An initial expecta- 
tion from the Liouville pathway suggests that this term 
should lead to a cross peak. Indeed, the expression 



I(a,T,C) = e[a]6[C] 
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(18) 

has different oscillation frequencies during a and (. The 
2D spectrum is displayed in Fig. [12] this term does give 
rise to a cross peak under nonrephasing conditions just as 
expected from the Liouville pathway and from femtosec- 
ond 2D ES. Similar to most of the other terms in this 
topological class, there is no dependence on time-interval 
r, and again the excited-state lifetime is related to the 
peak amplitude. 

The fifth and final term we consider in this class is 
Gmpu- This is one of the most unusual unrestricted 
terms. Expressed in the variables of the experiment, this 
term becomes 



I(a,T,C) = e[-T}e[C + a-T] 



PoM 2 |M2l 2 /. 2 
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The Heaviside functions for this term need to be consid- 
ered carefully. This term is nonzero when n > 0, meaning 
field Ea before field Elo, and when r < 0, meaning field 
Eb before field Eq- In other words, there is not a strict 
requirement on the delays internal to the pump and probe 
sequences as long as the pump sequence 'happens' before 
the probe sequence. 
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FIG. 13. The result for term Gmpu- (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (bottom) Absolute value I^ 2D ES 
spectrum. This term oscillates during r at the difference fre- 
quency of the two single-exciton states due to polarization 
interference, and the oscillation decays at the difference of 
excitonic dephasing times, j eig — lge 2 - 



The 2D spectra for this term under four different tim- 
ing conditions are displayed in Fig. 



13 The Liouville 



pathway suggests that, under rephasing time-delay con- 
ditions, this term should give rise to a cross peak that 
should decay slowly with no oscillatory behavior; in I^ 4 ) 
2D ES this term leads to a diagonal peak with an os- 
cillating amplitude during time period r which dephases 
(very slowly) at a rate given by r ) ei g—lge 2 • Under certain 
dephasing conditions the signal could even grow in am- 
plitude as a function of r. As a reminder, here we have 
chosen j eig = r ) ge21 thus the oscillations neither grow nor 
dephase during r. 

At first glance, the oscillations during r appear to be 
the quantum beating observed in many femtosecond 2D 
ES measurements. The oscillations are not true quan- 
tum beats; they are a manifestation of polarization inter- 
ference, a topic considered in the femtosecond literature 
throughout the 1990s regarding how a pair of uncoupled 
transitions can lead to coherent oscillations in certain 
one-dimensional four wave-mixing measurements such as 
self-diffraction |158H162j . It may be worthwhile to revisit 
and investigate this distinction between the femtosecond 
and noisy-light methods in future work. 

There are 32 nonzero terms in this class, and their con- 
tributions can be categorized by the four types of Heav- 



iside restrictions, 0[±cr]0[£] and 0[±r]0[^]. In Fig. 14 



we show the total results under six different pulse-timing 
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x = -0.2ps r = Ops x = +0.2ps 




FIG. 14. The total result for the unrestricted topological class 
for three values of r. Each spectrum is the absolute value of 
the sum of all the terms which contribute under the spec- 
ified time ordering, (top) The nonrephasing contributions, 
(middle) The rephasing contributions flipped into the (+,+) 
quadrant, (bottom) The total I (4) 2D ES. 



conditions, £ from to 0.5 ps, ±<r from to 0.5 ps, for 
r G {0, ±0.2} ps, and we show the total spectrum (sum 
of rephasing and nonrephasing) for each of the three r 
values. We chose these values based on expected exper- 
imental conditions. Each spectrum is the absolute value 
of the sum — not the sum of the absolute value — of all the 
terms which contribute under the specified time ordering. 

Some general features for the unrestricted topology 
terms are as follows. The terms can lead to diagonal and 
cross peaks. A few terms lead to peaks that oscillate due 
to polarization interference during time-delay interval r, 
but most do not oscillate or even decay as a function of r. 
No true coherent quantum beats are observed for terms of 
this topology because color locking eliminated the terms 
which could have done so. There are no two- quantum 
pathways in this topological class that contribute to the 
I^ 4 ) 2D ES signal. 



B. Singly Restricted Examples 

In this topological class we study terms G//i2/3i, 
Gn22pi, Gu22aH, C//i2/3i? C//i2mi? an d T^vv p%- We be- 
gin with G ni 2(3i, 

i(*,T,Q = e[-*]e[c + *] (20) 

x _ 



The product of Heaviside functions, 6[— cr]Q[( + a], 
means that even though this term originates from a non- 
rephasing pathway, it is active only during the rephasing 
pulse time ordering, see Fig. [l5^a). The restriction on £ 
originates from its coupling with a in the second Heavi- 
side function, where a is restricted by the first Heaviside 
function. Although the maximum amplitude in the time 
domain is quite small, Fourier transformation results in a 
peak in the 2D spectrum in Fig. [l5^c) with a respectable 
amplitude because the oscillations persist for much longer 
than the optical dephasing time along £ = — a. The peak 
will gain amplitude as either or both of the variables are 
scanned to longer times. The peak has a very small anti- 
diagonal linewidth. Close inspection reveals that there 
are two peaks, one along the diagonal at each exciton 
frequency. This term has no r dependence. 

We next study term Gj/22/3*, which is very similar to 
the previous term except it was of the pathological type, 
which lead to the expression 

/(a, r, C) = e[-cr]6[C + cr}(a + £) e - <n «2n(*+O e +in M * 

(21) 

The resulting 2D spectrum is shown in Fig. [16] This 
is one of the highest-amplitude terms in this topologi- 
cal class. This term contributes only along the diagonal 
near the higher-energy exciton state and is independent 
of time-delay variable r. This peak will also gain ampli- 
tude as the delays a and ( are scanned to longer times. 

The third term we consider is Gu22aiii where elimina- 
tion of time-delay variable k leads to 

/(a, r, C) = e[-r]6[C + <r](C + a)e-^ (c+a) e +i ^ r . 

(22) 

Here the restraint on a is lifted, and it is replaced with the 
less-stringent restraint on r. This means that the signal 
persists throughout the surface of the a and £ variables, 
see Fig. pT^b). Like the previous terms, the amplitude is 
small but the signal can be acquired out to large values 
of C and a. Here we display the resulting rephasing and 
nonrephasing spectra in Fig. [l7^c) and (d), respectively. 
The shape of the peak in the nonrephasing spectrum is 
very similar — albeit at reduced amplitude — to many of 
the peaks that originate from terms in the unrestricted 
topological class. 

The next term we consider is the first nonzero term 
that originated from a C type response function. The 
time-domain data and the 2D spectrum for Cj/i2/3z> 

/(a, r, C) = B[-a]6[C + a)(a + ^-^.(^+0^0^^ 

(23) 

are shown in Fig. [l8^a) and (c), respectively. 

The quantum beats of interest are convolved into <r, 
which leads to a shift of the peak into the off-diagonal 
position rather than coherent oscillations as a function 
of time-delay variable r. The signal amplitude is quite 
small, and since this term does not persist at long scan 
times, it will remain at a reduced amplitude in the fre- 
quency domain. There is no dependence on time-delay 
variable r even though this a coherence pathway term. 
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FIG. 15. The result for term Gn\2$i- (a) The real part of 
the signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value 2D ES spectrum. This term con- 
tributes only along the diagonal at both exciton frequencies, 
and it is independent of time-delay variable r. 



The fifth term of note in this topological class is 
CiiyiolU, another C type term, where elimination of k 
leads to 

I(a, r, C) = e[-r]6[C + cr}(a + £) e -^ 15 (<x+C) e +^ eie2 r_ 

(24) 

This is one of the most interesting terms of the experi- 
ment since it has true quantum beats during time period 
r. This nonrephasing pathway, under time-delay condi- 
tions equivalent to a nonrephasing scan, appears almost 
identical to femtosecond measurement: it contributes to 
the diagonal peak in Fig. [l9^d) and it oscillates as a 
function of r. In the noisy-light measurement, however, 
this term also contributes under the standard rephasing 
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FIG. 16. The result for term Gn22^%- (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (bottom) Absolute value I^ 2D ES 
spectrum. This term contributes only along the diagonal near 
the higher-energy exciton state and is independent of time- 
delay variable r. 
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FIG. 17. The result for term Gii22&ii- (a) The double-sided 
Feynman diagram, the WMEL diagram, and the FTC dia- 
gram for this term, (b) Time-domain signal, (c) Absolute 
value I^ 2D ES nonrephasing spectrum, (d) Absolute value 
I^ 4 ) 2D ES rephasing spectrum. This term contributes only 
along the diagonal near the higher-energy exciton state and 
is independent of time-delay variable r. 
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FIG. 18. The result for term Cii2/3i- (a) The real part of the 
signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value 2D ES spectrum. This term 
contributes only along the diagonal near the higher-energy 
exciton state and is independent of time-delay variable r. 



time-delay conditions and when £ < 0, both of which are 
not expected from femtosecond measurements, see Fig. 
191)). We show the case of standard rephasing condi- 
tions in Fig. [l9jc). Unfortunately the contribution from 
this term is extremely weak; note the maximum ampli- 
tude in Fig. 19 ^d) is on the order of 10 -3 compared to 



many of the unrestricted terms that have amplitudes on 
the order of the 10 7 . 

The final term of note in this class is Tivvpu a term 
which originates from a two-quantum pathway, 



I(<T,T,() = e[<T}e[(-a](<T-()e 



(25) 
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FIG. 19. The result for term Ciii2aii- (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (bottom) Absolute value 2D ES 
spectrum. These peaks oscillate as a function of time-delay 
variable r. 



teresting in femtosecond 2D ES is convolved with a one- 
quantum frequency during a. This leads the peak to 
appear at the cross-peak position, but not at the two- 
quantum frequency. However, this term and the ones 
similar to it have an unusual, twisted peak shape, see 
Fig. 20 'c). Similar to so many others, this term does not 
depend on the value of r. 



There are 40 nonzero terms in this class, and their 
contributions can be categorized by the four types of 
Heaviside restrictions, 6[— cr]&[( + a], 9[a]6[C — a], 
9[-r]e[« + r], and 9[r]6[^ - rj. In Fig. [2l] we show 



The two-quantum frequency that makes this term in- 



the total results under the same six different pulse-timing 
conditions used in the unrestricted topological class. The 
rephasing terms are about an order of magnitude stronger 
than the nonrephasing terms and thus dominate the to- 
tal response. We have chosen to not present any results 
under which ( is negative. 

Generally, terms in this topological class are much 
weaker than terms in the unrestricted class, and they 
can lead to diagonal and cross peaks. A few terms lead 
to peaks that oscillate due to true quantum beats, but 
they are quite weak. Most terms do not oscillate or decay 
as a function of r. There are no rephasing pathways in 
this topological class that contribute to the I^ 4 ^ 2D ES 
signal. We have not shown the spectra that result from 
the signal that occurs at negative values of time-delay 
variable (. 
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FIG. 20. The result for term T 21 > 1 >p i . (a) The real part of 
the signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value 2D ES nonrephasing spectrum. 



C. Doubly Restricted Examples 

The first term we consider in the doubly restricted 
topological class is the coherence term Cinpi- Its ex- 
pression is 

1 = e[C]9[-a]p | Ml | 2 |^ 2 | 2 J 2 ( e[;j + c] 

x e +i^ e2ei cr e -iQ e2 g(c+cr) _ e +in gei a e -in e2g ( 
+&[-a - C] 

e -in e2ei C e +in gei (C+v) _ e +in gei a e -in e2g c | ^6) 

and the time-domain signal and rephasing 2D spectrum 
for this term are presented in Fig. |22|a) and (c), respec- 
tively. This term is only active under typical rephasing 
conditions when ( > and a < 0. The peak is located 
in the off-diagonal position, but as the expression shows, 
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FIG. 21. The total result for the singly restricted topological 
class for three values of r. (top) The nonrephasing contribu- 
tions, (middle) The rephasing contributions flipped into the 
(+,+) quadrant, (bottom) The total I (4) 2D ES. 



this term does not oscillate as a function of time-delay 
variable r. The shape of the peak is quite unusual due to 
the symmetry in the time-domain signal about the line 
£ = a. The amplitude of this term is very small, and 
increasing the size of the temporal scan range will not 
increase its amplitude. 

The expression for the second term of note in this class 
is C/i2cm, whose expression is 

7 = e[ K ]e[-T}p \^\WP i [r + k] 

h^(Q e2 g ^g 2ei -|- Qg ei ) 1 

+ift eoei T e -iQ e2g (K+T) _ e +in gei r e -in e2g K 



x e 



where, for brevity, in this case we have chosen not to elim- 
inate time-delay variable k from the expression. Never- 
theless, the simulation was performed in the space of <r, r, 
and (. The signal and resulting 2D spectrum are shown 



in Fig. 23 Although this term does contain quantum- 
beat oscillations in the factor of e _ ^ e i e 2 r in the first part 
of the expression, the signal is very weak. Its contribu- 
tion to the 2D spectrum in principle could be increased 
by scanning time variables a and ( to larger values. The 
resultant peak has a very small anti-diagonal line width 
reminiscent of several terms in the singly restricted topo- 
logical class. 
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FIG. 22. The result for term Cn2p%- (a) The real part of the 
signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value 2D ES spectrum. 




FIG. 23. The result for term Cn2 a u- (a) The real part of the 
signal in the time domain, (b) The double-sided Feynman 
diagram, the WMEL diagram, and the FTC diagram for this 
term, (c) Absolute value I^ 2D ES spectrum. 



We turn to term G 112/3% next 

f= e[C]9[-a]po|Mi| 2 |^| 2 Jo 2 f 
fc 3 (^e 2 g ~ ^gg + ^gei) 1 



9[a + C] 



x e 



+in gg (T -ifle 2 g(C+0-) _ +in gei (T -ifl e2 gC 



+0[-a - c] 



>. (28) 



This term is an example of one of the stronger terms in 
this class, primarily because the signal persists about the 
line = — a to large values of either variable. Moreover, 
the peak amplitude in the 2D spectrum is large because 



the signal is compressed about the diagonal. There is 
the a second portion of the signal — small and rapidly 
decaying — that leads to the weak cross peak. 
The final example is term T^'i'm, 



I= 9[c]e[g] Po |MilW/ 2 
h 3 (n eig -nfg + n e2g ) 



-iflfgcr -in e29 (C-cr) _ p -iQ eig (T -iQ e2 gC 



+Q[a - C] 



-itof9<p-itoe ig ((r-t) _ p -iQ ei9 cT -iQ e29 C 



>. (29) 



This term originates from a two-quantum Liouville path- 
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FIG. 24. The result for term Gn2/3i- (top) The double-sided 
Feynman diagram, the WMEL diagram, and the FTC di- 
agram for this term, (bottom) Absolute value 2D ES 
spectrum. 



way, and it leads to a small- amplitude cross peak in the 

As with the other similar 



2D spectrum, see Fig. 25 



terms, this term has no r dependence, it is very small in 
amplitude, and its unusual peak shape is caused by the 
peculiar response in the time domain. 

There are 40 nonzero terms in this class, and their 
contributions can be categorized again by the types of 



Heaviside restrictions. In Fig. [26] we show the total re- 
sults under the same six different pulse-timing conditions 
used in the previous topological classes. The response 
under rephasing conditions is about an order of mag- 
nitude stronger than response under nonrephasing con- 
ditions (where the two-quantum pathways contribute), 
and thus the rephasing contributions dominate the total 
doubly-restricted response. 

Generally, terms in this topological class are much 
weaker than terms in the unrestricted class, and they can 
lead to diagonal and cross peaks. Most terms, including 
all the large- amplitude rephasing contributions, appear 
as narrow stripes along the diagonal. A few terms lead to 
peaks that oscillate due to true quantum beats, but they 
are again weak. There are no one-quantum nonrephasing 
pathways in this topological class that contribute to the 
l( 4 ) 2D ES signal. 



D. Total Signal 

An experiment would measure the total signal — 
meaning the sum of all 128 terms — regardless of topol- 
ogy. Thus we could present the total signal expected for 
typical rephasing and nonrephasing experiments under 
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FIG. 25. The result for term T 12 , V(xi . (top) The double- 
sided Feynman diagram, the WMEL diagram, and the FTC 
diagram for this term, (bottom) Absolute value I^ 2D ES 
spectrum. 
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FIG. 26. The total result for the doubly restricted topo- 
logical class for three values of r. (top) The nonrephasing 
contributions, meaning the contributions to the total signal 
under temporal scan conditions analogous to a femtosecond 
nonrephasing 2D measurement, (middle) The rephasing con- 
tributions flipped into the (+, +) quadrant, (bottom) The 
total I (4) 2D ES. 
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erties of the light used to probe light-harvesting proteins 
in femtosecond measurements. 



E. Two-level System I (4) 2D ES 

In the process of solving the Bloch four-level system, 
we have also solved a simpler problem, I^ 4 ^ 2D ES for the 
two-level system. Many signal components for the Bloch 
four-level system do not exist for this system. In the 
end, only two G type and two P type Liouville pathways 
remain for a and f3 types. This leads to a total of only 
sixteen FTC diagrams. For brevity, we do not present 
the expected signal. 



V. CONCLUSIONS 

Linear absorption spectroscopy is, formally, a coher- 
ent measurement technique. A linear absorption spec- 
trum can be measured using a femtosecond laser and a 
suitable detector. Yet, it is well known that this mea- 
surement does not require a coherent light source, and 
in fact incoherent light sources are extremely common 
for measurements in the visible and infrared portions of 
the spectrum. This is a powerful reminder that the co- 
herence property of the light source plays a role less im- 
portant than perhaps expected in coherent spectroscopy 
measurements. 

Our results suggest that there is value in attempting 
to measure a 2D spectrum using incoherent light. To be 
sure there will be a number of challenges to collecting 
the I^ 4 ) 2D ES data, especially if the goal is to observe 
coherent quantum beats. First, the interpretation of the 
experiment is not as straightforward as the femtosecond 
case because of the cw nature of noisy light. Second, 
methods to suppress the signals arising from terms in the 
unrestricted topological class will be needed to resolve the 
terms that produce quantum-beat signals clearly. Third, 
a noisy-light source with a significantly broader spectrum 
than any used to date will be needed. 

We hope that the connections between I^ 4 ^ 2D ES and 
its femtosecond analogue will provide new insights into 
the interpretation of the corresponding femtosecond ex- 
periments. The results presented here show that I^ 4 ^ 2D 
ES should contain diagonal and cross peaks similar to 
femtosecond 2D ES, and they show that quantum beat- 
ing persists in the third-order nonlinear response even 
under incoherent excitation. This work addresses in an 
experimentally verifiable manner the important question 
of whether or not coherent quantum beating between ex- 
citon states is an artifact arising from the coherence prop- 



the various values of r used above. However, inspection 
of the relative magnitudes shows that the entire class of 
singly restricted terms and the entire class of doubly re- 
stricted terms amounts to total signal on the order of 10 4 , 
whereas the unrestricted amplitudes are on the order of 
10 8 . Thus the unrestricted terms completely dominate 
the response, and the total signal simply replicates those 



presented in Fig. 14 
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VI. APPENDIX 

The FTC diagram analysis presented in this work 
greatly reduces the amount of mathematical effort re- 
quired to produce an expression for each component of 
the final signal. Although there is a formal isomorphism 
between the mathematical expressions and the FTC di- 
agrams |137[ 138 j, it is informative to include selected 
analytic calculations in addition to the FTC diagram 
analysis. In this appendix we demonstrate the lengthy 
integration and algebra required for four of the 128 sig- 
nal components, one of which is zero. The approach is to 
insert the response function and electric fields into Eqn. 
[5] and the goal is to arrive at the expression in the Tables 
that correspond to the terms G/nim? Gnvifii-> Cn2aii, 
and Cn2ai- We reiterate that while perhaps new to read- 
ers more familiar with femtosecond-pulse measurements, 
the assumptions and protocols used here are based in al- 
most thirty years of experimental and theoretical work 
on noisy-light spectroscopy. 



A. Term Gim a i 

This is one of the simplest terms. The response func- 
tion for Gmiai — and in fact for all four Gnu terms — is 



In the second step we reduce the four-point time cor- 
relator into a pair of two-point time correlators via the 
assumption that the fields obey circular complex Gaus- 
sian statistics. This leads to two separate terms: G/nim 
and Gj-jiicm- We select the two-point correlator pair that 
corresponds to Gmim> 

I(ct,tX)= I dt 3 I" dt 2 [ 2 dt 1 {E A El){E c El ) 
J— 00 J — 00 J— 00 



' ih 3 



-iCl ei g(t2—tl) 



■iQ eig (t-t 3 ) 



(32) 



Notice how the upper limits of integration are not +00 
but rather the next field-matter interaction because — for 
the particular Liouville pathway of interest — the field- 
matter interactions must have occurred in a specific or- 
der. We then assign each field to a particular interac- 
tion time according to the time-ordered Liouville path- 
way depicted by the appropriate diagrams in Fig. [3] As a 
nonrephasing pathway with an a subscript, field 22^ (£1), 
Esfo), and Ec(ts). Inserting the expression for each 
field according to Eqns. [I] into Eqn. 32 leads to 

rt pt 3 rt2 

dts / dt 2 / dt\ 



x{E Q p{t l )e- iujtl E^{t2 - (7)e +iuJ ^ 2 - a) ) 
x(E p(t 3 +T- a)e- iu;(t3+r - a) £; V(s - K)e +iu ^-^) 



x 



PolMil 



m 3 

xe -in gg (t 3 -t 2 ) e - 



xe -i£i gg (t3-t 2 ) e -itte ig (t-t 3 ) ^go) Reorganization yields 



where the transition dipole fi\ corresponds to the label 
m Fig. |2) In noisy-light spectroscopy, in addition to the 
response function, we must account for the quasi-cw and 
stochastic natures of the fields. There are four major 
steps. We first insert the response function into Eqn. |5j 

/(<t,t,C) = / dt s I 3 dt 2 [ 2 dt^EAElEcElo) 
J— 00 J— 00 J— 00 

x PoK e _ ffieis (i 2 - tl ) 
ih 3 



ih 3 

X(j>(tl)p*(t2- 
Xe ~iu)t 1 e +iu{t 2 



-ti) 

19 (t-t 3 ) 



ft 3 

dts I dt 2 

-00 



(33) 



cr))(p{t 3 + T-a)p*(s-K)) 

— o-)—iu}(t3+T—o-)+iu}(s — K,) 



Xe -iQ ei g(t 2 -t 1 ) e -iQgg(t3-t 2 ) e -ifl ei9 (t-t 3 ) ^ 

where Iq = EqEq. In the fourth major step, we approx- 
imate the noisy fields as being fully incoherent — their 
coherence times goes to zero — which allows us to replace 
the two-point correlators with (5-functions, 



xe 



-iQgg(t 3 — t 2 ) ~' iVL ei g(t— 1 3 ) 



(31) 



A)M 4 / 2 



Notice how the signal, /, is a function of the experimen- 
tal time-delay variables <r, r, and £. The interaction-time 
variables internal to the response function, ti, t 2l £3, and 
t, will disappear under integration. In femtosecond 2D 
ES there is a clear correlation between the experimen- 
tal time-delay induced by, for example, a retroreflector, 
and the time period between field-matter interactions. In 
noisy-light, we must account for all possible interaction 
times for each experimental set of time-delay variables. 
While we must account for all possible interaction events, 
the time-delay variables are under experimental control. 



dU 



*3 



dU 



dt\ 



t 2 + a)S(ts + t — o- — s + k) 



Xe -iujt! e +iuj{t2-<r) e - 



■iuj(t 3 +r — a) e +iuj(s — k) 



xe 



ji^-t!) e -iQgg(t 3 -t 2 ) e -iQ ei g(t-t 3 ) ^ ^ ^ 



We are now almost in a position to carry out the integra- 
tion. However, one must carefully inspect the arguments 
of the (^-functions because blind analysis can lead to er- 
rors that are difficult to catch. Care must be taken to 
understand when the product of J-functions leads the in- 
tegrand to be zero. This subtlety is better exposed by 
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making a change of variables from absolute interaction 
times to interaction-time intervals, T n = t n+ \ — t n : 



h=t- T 3 , 

t 2 =t-T 3 -T 2 , 

t 1 =t-T 3 -T 2 -T 1 . 



(36) 
(37) 
(38) 



This, in addition to setting t = s in anticipation of the 
heterodyned I^ 4 ^ 2D ES experiment, results in 

„ I.. |4 t2 poo poo />oo 

xS(a - Tx)(J(C - T 3 )e-^( £ - T3 - T2 - Tl ) 

Xe +iuo(t-T 3 -T 2 -a) e -iuo(t-T 3 +T-a) e +iuo{t-K) 



xe 



-iQ eig T 1 —iQ gg T 2 —ifl eig T 3 



(39) 



The negative sign results from the change of integration 
limits. Simplification of the above expression yields 



-ih 3 

xSia-T^SiC-Ts) 



xe 



-iu;(T+rv-T 3 -Ti) 



, e -iQ eig T le -iQ gg T 2e -iQ eig T 3 ^ 



(40) 



What makes this term simple is that the three integrals 
are independent, allowing us to reorganize the expression 
to 



L 



dT 3 S(C-T 3 )e- i ^- UJ ^ 

/o 

pOO 

x / (lT 2 (~ iil!l!/l2 
Jo 



POO 

x / dT^cr-T^e-^^-^K (41) 
Jo 

We can work in any order to solve the integrals, and 
we choose to first solve the T 2 integral. The indefinite 
integral is 



/ 



dT 2 e 



-iQ gg T 2 = _J_ e -in gg T 2 ^ 
^99 



(42) 



Since £l n 



-iTgg, where T gg is a positive value, the 



function decays and the definite integral yields 

oo 



n. 



99 







ft 



(43) 



99 



We now handle the T\ and T3 integrals using the generic 
expression, 



/>oo 

/ dT x S(z-T x )e- l(Q -^ T *. 
Jo 



(44) 



Due to the (5-function, the integrand — and therefore the 
integral as well — is nonzero only when T x = z. Thus 



f 

j — c 



dT x 5(z - T x )e-^ n -^ Tx = e -W-«,)z^ ^ 



which is also true when the lower limit of integration is 
instead of —00 as long as z was a value inside the limits of 
integration: < z < 00. If instead z < 0, the integration 
would yield 0, not the given expression. In other words, 

POO 

/ dT x S(z - T x )e-^-^ Tx = e-^-^GIz], (46) 
Jo 

where 9 is the Heaviside step function. Using this result, 
the total expression for this term becomes 

Key r n- gokfjJ^ZL -^fr+o 
xe -i(fie ls - W )ce^] 

xe-^i'-^eM, (47) 

where we used the identity t+k = cr + (. We can simplify 
this expression to its final form 

7(<t,t,C) = e[a]6[C] P °lso 4 ^ e- in ^e- in ^ a . (48) 



B. Term Gmzpi 

Performing all initial steps, the change of variables, 
and setting t = s, the intensity for this term is 

/( g , T , C ) =e -M^C) Pol^ 2 K| 2 / 2 



/>00 /'OC /'OC 

< / dT 3 / cff 2 / dTi 
Jo Jo Jo 

c(J(-r 3 -r 2 + c-ri)(j(r2 + (7) 

, e +iw(Ti+T 3 ) 

, e -in eig T 1 e -iVL gg T 2 e -in e29 T 3 ^ 



(49) 



This term is slightly more complicated than the first term 
because one of the (5-functions contains multiple integra- 
tion variables. This means we need to consider the inte- 
grals in consecutive order. Using the generic expression 
above, the T\ integral yields 

pOO 

\ dT x b (-T 3 - T 2 - Ti + C) e - i{n ^-" )Tl 
Jo 

= e -^(^ ei9 -u;)(-T3-T 2 +C)e[_ T3 _ T2 + Q (50) 

This leads to a T 2 integral which is slightly different than 
the previous terms. However, because of the second 5- 
function, this term is also straightforward 

pOO 

/ dT 2 5 (a + T 2 ) e -i(f2 99 +"-f2e lfl )r 2@ j_ T3 _ r 2 + C] 
J 



e +i(n gg +u>-n eig )* Q ^_ T3 + (T + QQ[_ a ^ 



(51) 
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After these two integration steps, the total intensity is 
I(a,T,C) = 9He-^^(«.»-" w ^l^l 3j o 



poo 

/ dT 3 e[cr + C - T 3 ]e- i( ^^-^i^ T3 . 



(52) 

The T3 integral is unlike any we have encountered thus 
far due to the presence of the Heaviside function involving 
the integration variable. This Heaviside function has two 
consequences. First it forces the maximum integration 
limit of the T 3 integral to be ( + a, because if T 3 > a + £, 
then the argument of the Heaviside function would be 
negative and thus the result would be zero. Second, that 
change to the limit of integration requires that C + c'" > 0? 
which we can ensure by applying a Heaviside of a + £ to 
the resulting solution. Thus, 

/ dT 3 B[a + C - T 3 ]e- i( ^-^i^ T3 
Jo 

= 6[a + C] / dT 3 e- l ^- ne ^ )T3 , (53) 
Jo 



which can be solved to yield 

—i 



e[a + (] 



^e 2 g ^ei# 



The total intensity for the term is 

l = e [ _ g ] e -<n. 1 .C e -<(n. 1 ,-n,,) g goH!N!jo (55) 



X9[<7 + C] 



^e 2 # ^ei^r 



(1 - e -^ e29 -^is)(c+^ 5 



which simplifies to 



/(<x,T,c) = e[-<7]e[<7 + c] 



PolMll M ^0 



C. Term Cii2an 

Performing all initial steps, the change of variables, 
and setting t = s, the intensity for this term is 



/ — e -iw(r+«) 



+„) Po|/il| 2 |/i2| 2 / 2 



-ih 3 

POO poo 



POO POO POO 

/ dT 3 / <£T 2 / dr x 
Jo Jo Jo 



x^-T 3 -T 2 )£(t + 72 + 71) 

Xe +iu;(T 3 -T 1 ) 

Like the term above, even though the integration steps 
must be performed consecutively, the ^-functions make 



the first two integration steps straightforward. Integra- 
tion over Ti and T 2 in the usual manner yields 

T _ PoM 2 M 2 /q r +in aei r r -m e2 e 1 -n aei )K 
-ih 3 

POO 

x / ^T 3 e _i ^ e29_ ^ e 2 e i + ^ 9e i^ T3 
Jo 



xe^-T 3 ]e[-r-^ + r 3 ]. 



(58) 



At this stage we must again carefully consider how the 
Heaviside functions impact the range of integration and 
the total signal. We first look back at the original expres- 
sion and observe that the final signal will acquire O [k] due 
to the first (5-function, which requires k — T 3 — T 2 = 0, 
where T 3 and T 2 are positive values. Similar reasoning 
leads to the signal acquiring a factor of 0[— r]. 

We can make more progress by following the line of 
reasoning above for the first Heaviside function alone. 
It will act to change the upper limit of integration to k 
and to introduce the same factor of Q[k] we found above. 
Based on this analysis, the signal becomes 



/ = e[n]Q[-r\ pM ^f I K +i ^ 

—ih 6 



PK 

x / dT 3 e- i{ne ^- ne ^ n ^ )T3 e[T 3 -r - «]. (59) 
Jo 

Finally, we treat the second Heaviside function. We 
recognize that T 3 > k + r, but the lower limit of in- 
tegration depends on whether or not n + r is greater 
than or less than zero. Thus we must consider two cases: 
k + r > and k + r < 0. 

For the k + r > case, the lower limit of integration 
is what we initially expected, n + r, and the signal gains 
Q[k, + r]. For the n + r < case, the lower limit is 
zero (because the original value of the limit was T 3 = 0) 
and the signal gains Q[—(k + r)]. This yields a single 
expression that uses the Heaviside functions to treat the 
cases 

j = e[ft]e[_ r ] ^Ql^l 2 l^l 2 ^o 2 e+ ^ gei r e -i(^ e2ei -^ gei )^ 
ih 3 



6[r- 



+e[- 



J T + K 

K ] r 

Jo 



dT 3 e 



Integration over T 3 then yields 

i = eMe[-r] ^l 2 ^ 



h 3 (n 



e 2 g 



^e 2 ei + ^gei ) 



+in gei r 



xe 



-i(^ e2g -^ e2ei +^ gei )(«:+r) _ -i(^ e2g -^ e2ei +^ gei )^ 



+e[- 



l-e" 



(61) 
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This simplifies a bit to 

/= 9M9[-r]po|Mi| 2 |M2|% 2 



0[t + k] 



+iQe 2ei r —iQ e2 g(T+K,) _ +iQ gei r —iQ e2g K 



+Q[-T - k] 



(62) 



In one case the 6 [ft] Heaviside is redundant, and in the 
other case the 6[— r] Heaviside is redundant. For com- 
pleteness we include both in each case. The Heaviside in 
the first term ensures ensures k, > \r\ and corresponds to 
column 2 in Table IIIIl and the Heaviside in the second 



term k < \r\ and corresponds to column 3 in Table III 



D. Term Cii2 a i 

We turn now to term C/i2m> the partner FTC diagram 
to the previous. FTC diagram analysis predicts this term 
is zero due to color locking. The mathematical analysis 
is similar to the examples above but we choose to use a 
path that better exposes the effects of color locking. Like 
above, the intensity for this term is 

1=1 dt 3 I" dt 2 r dt^EAE^EcElo) 

J— oo J— oo J — oo 



„ poImi| 2 Im2| 2 „ 
x =iP e 



id gei (t 2 — ti) 



-iQ e2g (t — t3) 



(63) 



Now though, when we reduce the four-point time corre- 
lator into a pair of two-point time correlators, we select 
the two-point correlator pair that corresponds to the i 
subscript, which yields 



poM 2 M 2 / 2 

*<P(t2)p*(tl- 



dt: 



/*3 /**2 
dt 2 / dh 
-oo J — oo 



<r))<p(ta+T-a)p*(a-K)) 



xe 



-iQ gei (t 2 —ti) -iQe 2 e 1 (t 3 -t 2 ) -ifie 2 g(t-t 3 ) 



, (64) 



At this point it is convenient to deviate from the path 
taken for Cn2 a u and consider (p(t2)p*(t\ —cr)) which we 
can write as 



j I'OO rC 

wL dxi L 



(27T) 

{p(x 2 )p*(x 1 ))e 



dx 2 



-ix 2 t 2 ix\{tx-a) 



,(65) 



where p(x) is the Fourier Transform of p(t). Now invok- 
ing the Weiner-Khintchine theorem, Eqn. [7| we have 



(P(t2)p\h-C7)) 



dx\ 



dxo 



T(x 1 )S(x 1 - x 2 )e- iX2t2 e iXlitl - a \ (66) 
where T(x) is the spectral density of the noisy light. The 
5{x\— x 2 ) is the mathematical manifestation of color lock- 
ing as it represents the x—x correlated pair. Inspection of 
the Liouville pathway — the WMEL or double-sided Feyn- 
man diagrams — for this 'C term shows that x\ ^ x 2 . 
Consequently, 



<p(t 2 )p*(t 1 -a))=0. 



Thus 1 = and this term vanishes. 



(67) 
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